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Abstract 

We present the exact ground-state wave function and energy of the generahzed Hub- 
bard model, subjected to the condition that the number of double occupied sites is 
conserved, for a wide, physically relevant range of parameters. For one hole and one 
double occupied site the existence of the ferromagnetic ground-state is proved which 
allow one to determine the critical value of the on-site repulsion corresponding to the 
point of metal-insulator transition. For the one dimensional model the exact solution 
for special values of the parameters is obtained. 
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1. Introduction 

The Hubbard model is tlie generic model to describe correlations in narrow-band 
systems The on-site repulsion is due to the matrix elements of the Coulomb inter- 
action corresponding to the on-site Wannier states while the other matrix elements are 
neglected. However the on-site repulsion can be sufficiently strong and the values of the 
matrix elements corresponding to the pair of nearest-neighbor sites can be comparable 
with the value of the simplest nearest-neighbor hopping amplitude. The Hamiltonian 
of the model which is often referred to as a generalized Hubbard model contains the 
interaction terms of the fourth order in the electron creation and annihilation operators 
corresponding to the nearest-neighbor sites. The Hamiltonian is 

<ij><J <ij>rJ 

+ V riiUj + niin2i, (1) 

<ij> i 

where a = 1, 2 - is the projection of spin, rii^ = cf^Ci„, rii = riu + and < ij > 
denotes the pair of nearest neighbor sites. The generalized model (|l|) has been studied 
previously by several authors p, 0, ^ . The other models with a similar kind of hopping 
term where considered in ref's 0. Even for X,V <^ U the presence of interactions 
which directly couple nearest-neighbor sites should lead to new effects. For example the 
correlated hopping terms are believed to play an essential role in the formation of high 
-Tc superconductivity 0. In the present letter we consider the Hamiltonian (|l|) at X = t. 
At this particular value the model is much more simple than the conventional Hubbard 
model. Note that the corresponding region in space of parameters is quite realistic in 
view of the estimates of these parameters for different systems (for example, see Q 
and references therein). At the half filling the ground state of the model a.t X = t can 
be found exactly in any dimensions in a wide range of the parameters t, U, V. We also 
study the metal-insulator transition at the critical value Uc which can be found exactly 
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in our case in contrast to the Hubbard model where the Mott picture of metal-insulator 
transition is not directly applicable at least for a simple square or cubic lattice where 
presumably (in spite of the predictions based on the Gutzwiller approximation [Q) the 
system is the antiferromagnetic insulator at arbitrary U. Recently the same model was 
studied by Strack and VoUhardt with the help of supersymmetric representation. We 
show that the wave functions proposed are in fact the ground states of the model in a 
range of parameters which is much wider than the region found in ref . . 

First, we show that at the half filling the exact ground state of the model ai X = t 
and U > zmax(2t, V) {z is the coordination number of the lattice) is a highly degenerate 
state without the double occupied sites and the system is paramagnetic insulator. For 
the on-site repulsion U < Uc, where Uc is the critical value, which is not necessarily 
coincide with the obtained bound, the creation of holes and double occupied sites is 
energetically favorable. At this point the transition to the metallic state take place. In 
complete analogy with the Nagaoka theorem for the infinite-U Hubbard model [^] we 
prove that in the sector of Hilbert space with one hole and one double occupied site 
the state with the lowest energy is ferromagnetic. We show that at V < 2t the critical 
value is determined by the ground -state energy of the problem with one hole and one 
double occupied site and Uc = 2zt. At present time the stability of the Nagaoka state 
at finite concentration of holes is not proved. However it is supposed to be the correct 
ground state at sufficiently small concentration of holes |T^. Assuming the stability 
of the ferromagnetic state at finite concentration of holes and double occupied sites 
at low concentration of holes we find the density of holes below the point of metal- 
insulator transition. For a bipartite lattice we find another region of the parameters 
where the determination of the ground state is possible. For a square or cubic lattice 
ai U < 2zV — zmax{2t,V) the ground state is given by the state with the electrons 
occupying only one of the sublattices. In particular, at V > 2t the ground state is 
known exactly at arbitrary U. The transition between two different ground states take 
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place at f/ = zV. Finally we consider the generalized Hubbard model at X = t in one 
dimension. The dependence of the energy on the total spin for an arbitrary number of 
holes is studied. These results may be useful in the context of study of the stability of 
the ferromagnetic state for the infinite U Hubbard model in higher dimensions. It is 
shown that the model is exactly solvable at \^ = 0. We also comment on the behavior of 
different generalizations of the model (|]) which include the antiferromagnetic coupling. 



Brief description of the results presented in section 2 was given in ref. |11 



2. Generalized Habbard model in arbitrary dimensions 

Consider the model (|lD at the half-filling (n = ni + n2 = 1) at X = t. The Hamilto- 
nian 

H = -t J2 4aCjai^ - n-i^a - ^j-a) + h.C. (2) 
<ij>u 

+V riiUj + U"Yniin2i, 

<ij> i 

where h.c. stands for hermitean conjugate, conserves the number of double occupied 
sites N = J2i niin2i: [H, N] = 0. The eigenstates of H corresponds to a definite number 
of double occupied sites N which at the half filling coincide with the number of holes. 
First, let us prove that at f/ > zmax{2t,V) the ground state corresponds to N = 0. 
It is convenient to express the Hamiltonian in terms of the fermionic creation and 
annihilation operators of the holes {cf, q) and the double occupied sites {df, di) defined 
starting from the ferromagnetic state \F >= J^i^til^ >• The up-spin electrons are 
described by the Holstein-Primakoff hard-core bose operators {bf,bi). To obtain the 
interaction term ~ V one can make the following substitution: 

■riii = ni + Hdi, n2i = l-Ui- Ud- 

The Hamiltonian is 

H = -tJ2 i^Cj + dtdj){bpi + (1 - ni){l - rij)) + h.c. (3) 

<ij> 



+y i^ci - ndi){ncj - Udj) + UN + zVL/2, 

<ij> 

where Ud = cfci, = dfdi, rii = bfbi, and the constraint rid + + < 1 which is 
equivalent to the infinite on-site repulsion between the particles is implied. N = J2i ^di 
is the number of double occupied sites and the energy NU is due to the last term in 
Eq.®. 

a) Variational theorem. The upper bound for the ground-state energy Eq < 
zVL/2 {L is the number of lattice sites) of the Hamiltonian can be obtained using 
the variational wave function with singly occupied sites 

l'^>= 114 n 4,10 >, (4) 

where £ and C are arbitrary disjoint sets of lattice sites which together build up the 
total lattice. In the representation (|^) |0 > is the state without the fermions. It is easy 
to derive the lower bound for the energy Eq. The Hamiltonian (j^) can be considered 
as a matrix H^is where the indices a,f3 are enumerate the set of possible configurations 
of particles on the lattice a = {ii...iN\ji---jN\h---hi) where are the coordinates 
of holes, double occupied sites and hard-core bosons respectively (see Appendix). One 
can see that due to the Fermi statistics of c- and d- particles the non-diagonal matrix 
elements of corresponding to the kinetic energy term of the Hamiltonian are 
equal to ±t. Clearly, for the Bose statistics these matrix elements would be equal to 
— t. The diagonal matrix elements are determined by the second term of Eq.@. The 
following theorem can be easily proved. For each of the eigenvalues E of any hermitean 
matrix Ha/s at least one of the inequalities 

\E — Haa\ < X! l^aisl (5) 

is satisfied. In particular, for the restriction on Eq one should take the minimal value 
of Haa {—zVN). The right-hand side of this inequality is determined by the number 
of hopping processes allowed for a given state a. Since the total number of c- and d 
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- particles is 2N the maximal value of the right-hand side of Eq.(|^) is 2ztN. In this 
way we find the following lower bound for the ground-state energy: Eo{N) > {—2zt — 
zV + U)N + zVL/2. One can further improve this estimate in the following way. For 
example, the hopping of c- particle to the nearest neighbor site occupied by d- particle 
is not possible and one should not take into account these terms in the right-hand side 
of Eq.(P). However in this case the contribution —V to the diagonal matrix element 
Haa does exist. In the opposite case of isolated c and d- particles the binding energy 
is absent while the hopping processes are possible. Thus we obtain the following lower 
bound for the energy: 

^o(A^) > {-z max(2t, V) + U)N + zVL/2. (6) 

Taking into account the upper bound Eq < zVL/2 and Eq.(|D we see that at U > 
zmax{2t, V) the ground-state wave function does not contain the holes and the double 
occupied sites, = 0, and the energy is exactly equal to Eq = zVL/2. The ground state 
is 2^-fold degenerate and is given by Eq.(^). For comparison the bound for U found in 
ref. is 4zt + zV. The authors use the following method to obtain the lower bound for 
the energy. With the help of the operators 

the Hamiltonian can be represented in the form 

<ij>a i 

+ {U- Azt) niin2i + V Y ^i'^r (7) 

i <ij> 

Since the average of the first term in Eq.(0) over the ground state is positive-definite 
the bound U > 4zt + zV is obtained. Clearly for an arbitrary lattice, at t,V U and 
a small deviation t — X the degeneracy is absent and the model (|I|) reduces to an 
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effective Heisenberg model with antiferromagnetic coupling constant 4(t — X)"^ /U . For 
a bipartite lattice one has the antiferromagnetically ordered ground state. 

For U < Uc the creation of holes and double occupied sites take place. We shall see 
that aX V <2t the critical value Uc is determined by the value of the lowest energy for 
the state with one empty and one double occupied site (A^ = 1). This problem can be 
solved exactly in two and three dimensions. The Hamiltonian (Q) has the same form as 
an analogous Hamiltonian for the infinite-f/ Hubbard model away from the half filling. 
The only difference is the absence of d- fermions in the latter case. The hopping term for 
the bosons is absent as well and their interaction with the fermions has the same form. 
According to the Nagaoka theorem [§, for one hole the ground state of the infinite-?7 
Hubbard model is ferromagnetic (the total spin is maximal). In our case it is possible 
to prove that at = 1 i.e. for one hole and one double occupied site the ground-state 
is ferromagnetic. Following the original Nagaoka proof consider the hopping process in 
which the hole and the double occupied site starts from a given positions and come back 
to the same positions after a number of steps. The configuration of spins (the bosons) 
can be different in the initial and the final states. For one hole the corresponding 
energy-dependent self-energy part introduced in ref.[^ is positive. The same would be 
true for arbitrary number of holes if the holes would obey the Bose statistics in a sense 
of the Hamiltonian (|^). Clearly in our case the positivity condition is satisfied since for 
A^ = 1 the statistics of c- and d- particles which represent the two different species of 
fermions is not important. The other steps of the Nagaoka proof can be applied without 
modification (the presence of the attraction ~ V in Eq.(^ is either not important). 
Thus at A^ = 1 the state with S = Smax — 1 which is the maximal total spin for a state 
with one hole is the ground state of the system. In the representation (^) that means 
that the number of bosons is equal to zero for the projection of spin = —{Smax — !)• 
The ground-state energy of two particles has the form Eq{N = 1) = —2zt + U + a/L, 
where the last term is due to the interaction of c- and d- particles. Let us suppose that 
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the interaction part of the energy is positive a > (the corresponding condition for the 
parameter V will be found later). Then the critical value of U is 

= 2zt. (8) 

At U > Uc the ground state is given by Eq.@ and at U < Uc the creation of holes is 
energetically favorable. The critical value @ is the point of metal-insulator transition. 

Before solving the two-particle problem let us comment on the behavior of the system 
at finite density of c- and d- fermions. At present time the stability of the Nagaoka state 
for the infinite- repulsion Hubbard model at finite concentration of holes is not proved. 
It is supposed that the ground state is ferromagnetic at sufficiently small concentration 
of holes (for discussion see ref.|jlO|). In our case we will also assume the stability of the 
ferromagnetic state (5 = Smax — N) at sufficiently small density of holes and double 
occupied sites p = N/L. Then at U < 2zt and \2zt — f/| the density p — and the 
analog of the Nagaoka state is realized. The system of two species of interacting fermions 
at equal density p should be considered. At low density the ground-state energy can be 
evaluated as a series in the small parameter (|lnp|^^ and p^^^ respectively in two and 
three dimensions). At low density the dependence of the energy on p has the form 

E/L = -{2zt - U)p + ^o(p) + \a{p)p\ (9) 

where the second term is due to the Fermi statistics. In the lowest order £q{p) = intp'^ 
in two dimensions and £o{p) = (|6^/^7r^/^)tp^/^ in three dimensions. In the lowest order 
in p the function a(p) is determined by the two-particle scattering amplitude at low 
energy. The scattering amplitude can be expressed through the interaction energy of 
two particles of different species a/L in the finite volume L. In fact the low density 
limit of a(p) coincide with the value of the parameter a. The behavior of the function 
a(p) is different in two and three dimensions. In three dimensions it is the constant 



in the lowest order in p while in two dimensions a(p) = 87rt/| lnp| + 0(t/(lnp) ) [12 . 



The ground-state density of holes is determined from the condition of minimum of the 
expression (^). In the hmit U 2zt the minimum is determined by the first two terms 
in Eq.(^. For example, in two dimensions 

2zt-U ,^ , 

Po = (10) 

with the accuracy up to the terms of order ~ 1/lnpo- One can also take into account 
the leading corrections in po in this formula, po is small at small deviation of U from 
the critical value (^) in agreement with our assumption about the stability of the fer- 
romagnetic state. Although the description of the system (|^) at arbitrary number of 
bosons is not possible it will be shown that at Izt — U ^ the minimum of the energy 
as a function of p found from Eq. (^ is the correct ground-state energy of the system in 
agreement with Eq.(P). 

b) Solution of the two-particle problem. Let us proceed with the solution 
of the two-particle problem. The interaction potential (Q) contains the infinite on-site 
repulsion f/ — oo and the attraction of strength V at nearest neighbor sites. Let us 
consider the case of three dimensional cubic lattice. The ground state corresponds to 
the total momentum of two particles equal to zero. One can seek for the wave function 
in the form 

\^|J>=Y.Fik)ctdt,\0>, (11) 

k 

where = L^^^"^ Y.i^^^'^cf . The function (pH]) is the eigenvector of eigenvalue E if the 
function F{k) satisfies the Shrodinger equation 

(E - 2te,)F{p) = W{U- Vek-p) F{k), (12) 

k 

where ek = — 2 X]i=i cos /c^. Let us define the function J{k) = {E — 2tek)F{k). The 
equation (|l^) takes the form 



In order to solve the equation ([T3| ) for the ground-state energy Eq = —2zt + a/L let us 
extract from the sum the term with k = 0, which is most divergent for L ^ oo. We can 
also substitute the value 2teo (cq = —z) for E in the sum over /c 7^ in Eq.(|TB|) since 
in three dimensions the energy difference — eo is at least of order ~ L""^^^ and the 
interaction correction is of order ~ 1/L. We get 

Defining the new function r(A;) = aJ{k)/ J{0) we obtain from (JT^) the equation 

2tT{p) = U-Ve, + jY: ^-^''-^ m, (15) 
^ fc^o ^0 ^ 

which is nothing else but the equation for the scattering amplitude at zero energy r(0) = 

a. The sum can be replaced by the integral in Eq. (p!5D . The solution of the equation 

(|T5D can be represented in the form T{k) = Tq + efcFi. Substituting this function into 

Eq.(p!5|) we get two equations for the unknown constants FojEi. In the limit U ^ 00 

the result of the calculations for the simple cubic lattice {z = 6) has the form 

z{2t - V) 
^ ~ Wz-V{Wz -l)/2V 

where W = 0.2527 stands for the Watson integral 

2 rTT rTT rTT ^ 

W = III dkrdk„dk. 



(16) 



2{2txY J J -IT J -TV ^ ^ ^3 — COS k^ — COS ky — COS kz 
The expression (|l^) is valid at V < 2t. At V = 2t the amplitude a vanishes which 
indicate the existence of the two-particle bound state at \^ > 2t. Vanishing of the 
interaction correction to the energy aX V = 2t can be seen from analogy with the two- 
magnon problem in the ferromagnet where in the ground state the total spin should 
be maximal. The same conclusion can be made for two dimensions. To calculate the 
parameter a one should substitute the sum 

— ~ — In L 
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for W in the formula (pIBI). The corrections to the equation (|T^) are of order (1/lnL)^ 
and the terms of that order in Eq.(0) cannot be fixed. However the expansion in 1/ InL 
breaks down only at 2t — V ~ 1/ In L and at these values of V the perturbation theory in 
2t — V can be used since aXV = 2t the ground-state wave function is known: PgCq^c^o^IO > 
(the analog of 5* = Smax state in the ferromagnet; Pg is the Gutzwiller projector). As 
in 3D the interaction correction changes sign aX V = 2t and the bound-state solution 
appears at > 2t. 

c) Formal proof of the equation Finally it is necessary to show that, whether 
or not the Nagaoka state is realized at a given density p, the correction to the energy, 

E/L = -{2zt-U)p + £{p), (17) 

is strictly positive £{p) > Q and does not vanish in the thermodynamic limit at p 7^ 0. 
That means that the density is really small in the neighborhood of the point of metal- 
insulator transition. We have to obtain the lower bound for £{p) (p^. Let us modify 
the Hamiltonian (^ in such a way in order to decrease the corresponding ground-state 
energy. First, let us replace in Eq.(^) the repulsion V I](^ci^cj + ^dj'^di) at nearest neigh- 
bor sites by the attraction of the same form {V —>■ —V) and then make the substitution 
V ^ 2t so that the resulting interaction takes the form 

-2t {rici + ndi){ncj + Udj) + UN. 

<ij> 

Second, instead of fermions, consider the particles c, d, obeying the Bose statistics. Since 
for bosons the Nagaoka state is the ground state at arbitrary density and the Hamiltonian 
is symmetric in respect to the replacement c ^ d, the ground-state wave function 
(f){ii...iM\ji---jN) which is the totally symmetric function of its arguments coincide with 
the ground- state wave function of the Heisenberg ferromagnet in the representation of 
the Holstein-Primakoff bosons with the projection of spin = Smax —N and S = Smax- 
In fact, the corresponding wave function, which is the positive-definite and the totally 
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symmetric function, is the eigenstate of the modified Hamiltonian. Actually it is the 
ground-state of the modified Hamiltonian since for the wave function which changes sign 
the substitution 

0(il..iAr|ji...jAr) ^ |0(il..iiv| jl-.-jTv) | 

would lower the energy and the positive-definite eigenstate is unique because of the 
orthogonality condition. Note that these considerations can be a basis of a simple proof 
of the Nagaoka theorem both in our case and in the case of the infinite-U Hubbard model 
[P]. Therefore the lower bound for the energy is {—2zt + U)N and we find that S{p) > 0. 
Clearly at finite density S{p) does not vanish in the thermodynamic limit. Thus it is 
proved that at 2zt ~U the value of p minimizing the energy po ^ 0. Consequently, 
a.t V < 2t the point of metal- insulator transition is indeed given by Eq.(^). 

According to the Mott picture for large coupling U the density of states exhibits two 
bands with the centers separated by U. In the absence of electron correlations the width 
of each Hubbard band is zt, and the gap between the bands is expected to vanish at 
U/2zt = 1. The Hubbard bands are usually obtained in the framework of a special single 
- particle Green function decoupling approximation scheme proposed by Hubbard IQ , the 
Hartree-Fock - type approximation which was exact for zero interaction energy or zero 
band-width (which amounts to a specific decomposition of a certain Green functions and 
has no justification; for example see ref. |]l3l for a critical discussion). However, strictly 
speaking, in the conventional Hubbard model the very notion of the Hubbard bands is 
justified only in the limit U —>■ oo, since the number of double occupied sites is conserved 
only in this limit. In our model the number of double occupied sites is conserved and 
the notion of the Hubbard bands has a precise meaning beyond the framework of any 
approximation. Note also that our results can be used to explain the results of the 
numerical calculations for a small lattices [H . 

d) Ground state for a bipartite lattice. For an arbitrary bipartite lattice one 
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can find another region of the parameters where the determination of exact ground state 
is possible. For example consider the simple square or cubic lattice. The wave function 
which corresponds to a charge- density wave with maximal order parameter, 

lx>= 114410 >, (18) 

ieA 

where A is one of the sublattices, is an eigenfunction of H. The wave function (|18D can 
be used to obtain an upper bound for the energy: Eq < UL/2. To obtain the lower 
bound it is convenient to define the operators cu = (— l)*cij, C2i = (— l)*c^j. At the half 
filling the particle number corresponding to the new operators is the same. In terms of 
the operators Qo-, the state |x > is an antiferromagnetically ordered state with singly 
occupied sites. Since in terms of Qo-, the kinetic energy term has the same form as 
in Eq.(|^), in this representation the number of double occupied sites is conserved and 
the lower bound for the energy of the state with N holes can be obtained using the 
representation (|^) and the theorem (|^). The energy of isolated c and d- particles is zV 
and their interaction at nearest neighbor sites is —V. Thus the lower bound for the 
energy as a function of N is 

Eo{N) > {2zV - z max(2t, V) - U)N + UL/2. (19) 

One can see from Eq.([T9|) that the wave function ([18D is the ground state at f/ < 
2zV — zmax(2t, V). This ground state is unique apart from a twofold degeneracy due 
to the two sublattices and describes the nonmagnetic insulator. At V > 2t we obtain 
the condition U < zV. Since it was shown that a.t U > zV the ground state is given 
by Eq.(^), for a bipartite lattice at V > 2t the ground state is found for arbitrary U. In 
the Hubbard model (X = 0, V = 0) the function (^ is the ground-state wave function 
only for U = oo. In the generalized Hubbard model with X = t this ground state is 
already realized at finite U with U > zma.x{2t,V). Similarly, while in the extended 
Hubbard model (X = 0, V 7^ 0) the function (|T^ is the ground-state wave function only 
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for y = oo, in the generalized Hubbard model with X = t this ground state is already 
realized at finite V. 

3. Generalized Hubbard model in One dimension. 

Let us consider the generalized Hubbard model (0) in one dimension. Before consid- 
ering the model (0) let us study the ground state multiplicity (the value of the total spin) 
in the cases when the model is exactly solvable in one dimension. Namely we consider 
the model at V = and at 7^ in the sector with no double occupied sites. 

For the one dimensional system the two different cases can be considered: 1) the 
case of an open chain; 2) the case of the closed chain of finite length. Both cases are 
equivalent in the thermodynamic limit {L 00). First, consider an open chain or 
equivalently the chain of an infinite length. In this case the Nagaoka theorem is valid 
at arbitrary number of holes and double occupied sites since the Hamiltonian (0) is 
invariant under the transformation which changes the statistics of holes. In the other 
words for an open chain the holes can be considered as a bosons. That does not mean 
that the eigenstates with S < Smax cannot be degenerate in the energy with the state 
S = Smax- In fact all the eigenstates including the ground state are degenerate in the 
total spin S. For example let us consider the eigenstates of the Hamiltonian (|^) without 
the double occupied sites. At = in the case when both the numbers of c- and d- 
particles are not equal to zero we get the additional degeneracy of the eigenstates due 
to the two species of particles. Let us seek for the ground state wave function in the 
following form 

ij{ii, . .AnVi, . . . hi) = i^o{ii, ■ . .i7v)0(Ai, . . . \m) (20) 

where ia are the coordinates of c- particles {N = Nc) and A/3 are the coordinates of the 
spin bosons on a "supperlattice" which consists of Li = L — N lattice sites which are 
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not occupied by the holes {Li = is the number of electrons) 

N 

K = L-'^OiL-i/s), a = l,...M. 

(3=1 

If ipQ is the eigenstate of the Hamiltonian in the sector S = S^ax then the wave func- 
tion (^) is the eigenfunction of the Hamiltonian for an arbitrary function 0(Ai, . . . Am)- 
Thus the ground state is degenerate in the total spin S. Let us turn to the case 2) and 
see how this degeneracy is resolved at finite L. Note that unlike the open chain for the 
closed chain the spectrum depends on the statistics of particles. The coordinates on 
a superlattice can be defined in the same way by fixing the initial and the final sites 
of the chain {l^ = I, . . . L and = 1, . . . Li). After the substitution (^) we consider 
the functions ipQ and 0o extended to the infinite chain and subjected to the appropri- 
ate boundary conditions. In order to satisfy the periodic boundary conditions for the 
function ip ( pOD the functions i^oi^ii, ■ ■ - iN) and 0(Ai, ... Am) should satisfy the following 
boundary conditions: 

0(Ai,...A„ + Li,...Am) =0(Ai,...A„,...Am) a = l,...M, (21) 

and 

ipoiii, . . + L,. . An) = exp(ig)?/'o(H, i^v), (22) 

where the boundary conditions for the function ipQ are determined by the total momen- 
tum q corresponding to the function (p : 

0(Ai + 1, . . . Am + 1) = e'V(Ai, • • • Am). (23) 

The periodic function cf) is the symmetric function of its arguments which vanishes at 
Aq = A/3. At arbitrary V the function ipQ is determined by the set of the momenta /c^, 
a = 1, . . . A^. The energy is 

N 

E = —2t ^ cos ka- 

a=l 
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For instance at V = we have the free fermion determinant 



ipo = det [exp{ikaii3)] , ka = 27ma/L, 

where are the integers. The periodic function can be characterized by the set of 
the momenta Qa = 2nma/Li where rria are integers or half integers (see below). An 
arbitrary number of zeros = is possible which corresponds to the value of the 
difference S — S^. The total momentum q = J2a 2nma/Li. Thus we obtain the equation 
for the momenta ka'- 

= (24) 

As a basis in the space of the symmetric functions (0) one can choose the eigenstates 
of the Hamiltonian of the free hard core bosons (XX model) on a chain of the length 
Li. In this case rria are integers (half integers) for M- odd (even), and rria 7^ for 
a ^ f3. This automatically gives the eigenf unctions of the Hamiltonian (^. However in 
general this functions are not the eigenf unctions of the operator of the total spin S. In 
order to classify the eigenstates according to their spin one can choose the basis given 
by the eigenstates of the Heisenberg ferromagnet. In this case an arbitrary number of 
the momenta can be equal to zero and the non-zero momenta are determined by the 
system of the equations 



M 



^i,.L, ^ (_i)M-ig^p I 2z^^arctg(M« - n^j) | , = ^ctg(ga/2). (25) 



The total spin S = \Li/2 — M|, where M is the number of the non-zero momenta. 
The total momentum is g = 27rma/Li and we obtain the equation (p^ where 
are integers and rria are integers (half integers) for M- odd (even) (we assume L to be 
even). The same formulas could be obtained starting from the problem of Ng spinless 
fermions and the hard core bosons (upturned spins) on a " superlattice" consisting of 
the lattice sites occupied by the electrons. The only difference is the number of the 
momenta ka- a = 1, ---N^- In this representation the same results could be obtained by 
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taking the limit [/ oo in the exact Bethe anzatz solution of the ID Hubbard model 
[|T^ (for example see ref. |T^). In fact one can redefine the quantum numbers according 
to Ua ^ n^ — M/2, TTLa — > + Li/2 to obtain the equation with 71^- integers (half 
integers) for M- even (odd) and ma- integers (half integers) for (Li — M)- odd (even) 



in agreement with the results of ref . . 

From Eq.(^) the ground state energy as a function of the total spin can be found. As 
an example consider the splitting of the energy levels with S = Smax and S = Smax — 1 
for an even number of holes {L- is assumed to be even). Clearly the ground state 
corresponds to the values q = and q = n respectively and the value Eq = Eo^Smax ~ 1) 
is an absolute ground state energy. Thus we obtain the energy Eo{Smax) = Eq cos{it/L), 
where Eq = —2t (sin(7r/L))~^ sin(7rA^/L). For all 5* < Smax ~ 1 the minimal energy levels 
are nearly degenerate (i.e. the energy splitting is of the higher order in 1/L at large L) 
with the ground state energy Eq. The similar picture can be obtained for odd. In 
this case we find E{Smax) = Eq (M = 0, g = 0) in agreement with the Nagaoka theorem 
(A^ = 1). Clearly the same procedure (P(]|)-(P5|) could be performed at arbitrary V if 
the double occupied sites are absent or at V = and an arbitrary number of the double 
occupied sites. For instance in the first case the Bethe anzatz equations for the wave 
function ipQ ( pOD with twisted boundary conditions (^) should be used. 

Let us consider the model (0) at arbitrary U and arbitrary filling fraction n. It 
was proved that in the thermodynamic limit the eigenstates are degenerate. Thus it 
is sufficient to consider the eigenstates with S = S'' = Smax, which corresponds to 
the absence of the spin bosons. First, consider the model at = 0. At = the 
system is equivalent to the system of free fermions with an extra degeneracy due to the 
two different species of particles (c, d) with an infinite on-site repulsion. The density 
of double occupied sites p = N/ L is determined by the minimum of the energy of free 
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fermions with the total density no + 2p, 

2t 

Eo{p)/L = sin [7r(no + 2p)] + f/p, (26) 

TT 

where no = 1 — n is the concentration of holes in the limit of large U. Thus dX U > Uc 
where 

Uc = 4:t cos(7mo), 

the transition to the state with no double occupied sites, p = 0, takes place. Of course 
at no 7^ the point U = Uc is not related to the metal - insulator transition. Away 
from the half filling the analog of the ground state (|18|) i.e. the ground state without 
the single occupied sites is realized a.t U < —At. At V 7^ the model is not exactly 
solvable in the sector of Hilbert space with the double occupied sites. In the previous 
section it was shown that at arbitrary V < 2t and U > At we have the number of double 
occupied sites = (really N = at U > Uc where the critical value Uc < At). At 
these values the ground state is equivalent to the ground state of the Heisenberg chain 
with the anisotropy parameter V/2t in the sector with the projection of the total spin 
related to the number of holes. The spectrum of the anisotropic Heisenberg chain can 
be determined exactly with the help the Bethe anzatz. Thus in ID the ground state 
and the low energy excitations of the model (0) can be found exactly at U > Uc and an 
arbitrary n. At V/2t < 1 and V = 2t the spectrum of charge excitation is gapless and 
the system is a metal (if the number of holes no > 0). The same is true for V/2t > 1 
and no < 1/2. Note that away from the half filling the transition between the state with 
no double occupied sites and the state with no single occupied sites analogous to the 
states (§) and (0) does not take place at U = 2V under the condition V/2t > 1. If the 
concentration of holes is exactly no = 1/2 and V > 2t there is a gap in the spectrum of 



charge excitations [|T7|. Thus at the filling fraction ni = n2 = 1/4 the system undergoes 
another metal-insulator transition at the point V = 2t. In general the model (0) is not 
exactly solvable in the sector with the non- zero number of double occupied sites. 
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Although at the point of metal-insulator transition the ground state is ferromagnetic 
for two and three dimensional systems and degenerate in the total spin in one dimension 
the existence of the transition is not connected with the ferromagnetic order. In fact the 
metal - insulator transition is general phenomenon in the models with the kinetic energy 
term conserving the number of double occupied sites @. For example one can study 
the metal - insulator transition in the one dimensional models with an antiferromagnetic 



coupling which are exactly solvable in the absence of double occupied sites [p!8|] , |1T9| , [pO 
Apart from the term ~ X {X = t) and the on-site repulsion ~ U these models include 
the interaction of the form JJ2{SiSj — ^niUj) or JJ2{SiSj + ^niUj) at J = 2t. Although 
these models are not exactly solvable at 7^ at the half filling the existence of 
metal-insulator transition can be shown and the critical value of U can be found exactly 
Uc = 2zt In 2. The model which is exactly solvable at arbitrary N was proposed in ref. 
The Hamiltonian is the sum of the Hamiltonian of the tJ- model (modified to include 
the double occupied sites) and the permutation term of the form 

2t ^ {cfcjdj'di + Cj'Cidfdj). 

<ij> 

One can also change the sign J ^ — J to obtain the integrable model with the ferromag- 
netic ground state. The ground state of this model is equivalent to the ground state of 
the tJ- model [1^ with the up- and down- spin electrons replaced by c- and d- particles 
(that corresponds to the SU{2) "77- spin" symmetry of ref. [^). The concentration of d- 
particles should be found from the condition of minimum of the energy. 

4. Conclusion. 

In conclusion, for the model (^ it was shown that at the half filling and U > 
zma.x{2t,V) the ground state is given by Eq.(^). The problem with one hole and one 
double occupied site in the ferromagnetic background was solved. We proved that it is 
the lowest energy state in the sector of Hilbert space with one empty and one double 
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occupied sites. We established that at V < 2t the critical value of U corresponding 
to the point of metal-insulator transition Uc = 2zt. Under the assumption about the 
stability of the ferromagnetic state at finite concentration of holes the density of holes 
was found at V < 2t and U Uc- Finally, for a bipartite lattice with constant number 
of nearest neighbors at U < 2zV — z max(2t, V) the exact ground state wave function 
is given by Eg. (p!8D . A.t V > 2t the ground state is found exactly at arbitrary value of 
the parameter U . The transition between the states (|) and (|18D occurs at U = zV . At 
V <2t the exact ground state is found at U > 2zt (|) and U < 2zV - 2zt For 
the one dimensional model at V = Q the exact solution was presented. We have also 
studied the dependence of the energy on the total spin. Recently the part of the results 
obtained in Section 3 was independently obtained in ref . . 
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Appendix 

Following Nagaoka §, we introduce a set of orthogonal and normalized many-body 
wave functions which completely span the Hilbert space. We use the representation (|^) 
in terms of the operators of holes cf = C2i, double occupied sites df = cf^ and the 
overturned spins (hard core bosons) bj^ = Ci^C2i starting from the ferromagnetic state 
J2iC2i\0 >■ Because of the Fermi statistics we have to be cautious about the order of 
c- and d- operators in the definition of the states a, (3. Let for each site i of the lattice 
Ri to be an integer number Ri = 1, . . . L. We can define the following order among the 
lattice sites. Setting up the coordinate system we assign a pair of integer coordinates 
{ix, iy) to each site of the lattice. If for the two sites i, j ix < jx then we define Ri < Rj. 
When ix = jx the order is determined by their y coordinates: Ri < Rj for iy < jy. With 
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this order we now introduce the states 

l«>=<---4<.---<.^h---Clo>, 

where the order of the operators for a given set of the lattice sites a = {ii---iN\ji---jN\h---lM) 
is given by the condition 

Rii < ... < RiN, Rji < ... < RJn, Rh < ... < RIm- 

With this definition the non-diagonal matrix elements =< a\H\f3 > of the Hamil- 
tonian (0) are equal to —t when one hole (double occupied site) changes order with an 
even number of holes (double occupied sites). In the opposite case the matrix elements 
Hai3 are equal to +t. Clearly for the Bose statistics the matrix elements would be equal 
to —t. Since c- and d- particles are the two distinct species of fermions one can use 
the prescription given by the equation for \a >. For = 1 the non-diagonal matrix 
elements are also equal to —t which justifies our statement that the Fermi statistics of c- 
and d- particles is not important for the case of one hole and one double occupied site. 
Thus for the model (|^) the Nagaoka proof can be used without modifications for one 
hole and one double occupied site. 
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